R Exercises - G(r)

gofr <- function(d, dr = .01) {
    # Get a vector of all Euclidian distances
    dd <- as.vector(dist(d))
    # Make a histogram out of it
    dd.hist <- hist(dd,
        breaks = seq(0, max(dd) + dr, by = dr),
        plot = FALSE
    )
    # Get the r values
    r <- dd.hist$mids
    # Compute the normalization by the bin surface
    rlo <- r - dr / 2
    rup <- r + dr / 2
    nideal <- pi * (rup^2 - rlo^2)
    # Return the tibble containing r and G(r)
    tibble::tibble(
        r = r,
        gofr = dd.hist$counts / nideal
    )
}
# Example of use
ddr <- 1
d <- tibble(
    x = rep(seq(0, 10, ddr), 10 / ddr + 1),
    y = rep(seq(0, 10, ddr), each = 10 / ddr + 1)
)
ggplot(d, aes(x, y)) +
    geom_point()

G <- gofr(d, dr = .05)
library(plotly)
P <- ggplot(G, aes(x = r, y = gofr)) +
    geom_line()
ggplotly(P)